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Abstract 

We develop an analytical formalism to determine the statistical 
properties of a system consisting of an ensemble of vortices with ran- 
dom position in plane interacting with a turbulent field. We calcu- 
late the generating functional by path-integral methods. The function 
space is the statistical ensemble composed of two parts, the first one 
representing the vortices influenced by the turbulence and the second 
one the turbulent field scattered by the randomly placed vortices. 
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1 Introduction 



This work presents a study of the statistical properties of a system of vortices 
interacting with random waves. This is motivated by the necessity to describe 
quantitatively the statistics of a turbulent plasma in the regime where struc- 
tures are generated at random, persist for a certain time and are destroyed by 
perturbations. This regime is supposed to be reached at stationarity in drift 
wave turbulence when there is a competition between different space scales. 
These scales range from the radially extended eddies of the ion temperature 
gradient (ITG) driven modes to the intermediate a/ Lnps interval where the 
scalar nonlinearity is dominant over the vectorial one and with random fast 
decays to few scales where robust vortex-like structures are generated. 
In a physical model the description of these states requires to consider si- 
multaneously the long time response and the much faster events of trapping 
of the energy in small vortices, condensed from transient Kelvin-Helmholtz 
instabilities of the locally ordered sheared flow. 

The analytical model we propose here is essentially nonperturbative in 
the sense that the structures are explicitely represented and we take into ac- 
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count their particular analytical expression (or a reasonable approximation). 
It is in the spirit of the semiclassical methods, frequently used in field the- 
ory. In this sense this treatment presents substantial differences compared 
with the more usual approaches, aiming in general to obtain a renormaliza- 
tion of the linear response or to calculate the correlations of the fluctuating 
field using a closure method. When structures are present in the plasma 
any perturbative approach is inefficient, because it can only go not too far 
beyond Gaussian statistics, i.e. calculate few cumulants beyond the second 
one. The problem consists in that they take as the base-point for developing 
perturbative expansions the state of equilibrium of the plasma or a Gaussian 
ensemble of waves. Or, from this point it is impossible to reach the strongly 
correlated states of coherent structures. The idea of semiclassical methods is 
precisely to place the equilibrium state on the structures and then to explore 
a neighborhood in the space of system's configurations, to include random 
waves. This is the approach we will adopt in the present work. To have a 
tractable problem we assume the simplest case where the vortices are not 
created or destroyed dynamically. The final outcome from this approach is 
a list in which a particular dependence of a contribution to the correlation, 
e.g. the exponent fi in ((/x^)^ ~ k~'^, is associated to a particular contribu- 
tion in the physical model: statistics of the gas of vortices, interaction energy, 
nonlinearity of the physical model, etc. 

1.1 The formulation of the problem 

In agreement with the experimental observation and numerical simulations 
P] it has been found theoretically that the nonlinear differential equations 
describing strongly nonlinear electrostatic drift waves in two-dimensional 
plasma can have (1) turbulent solutions, consisting of very irregular fluc- 
tuations which can only be described by statistical quantities (irreducible 
correlations = cumulants); and (2) solutions that are cuasi-coherent struc- 
tures of vortex type, which are remarkably robust and for which we have in 
certain cases explicit analytical expression. Although there is no conceptual 
difference between these two aspects, we can say that we have two mani- 
festations of the nonlinearity: one is nonlinear mode coupling and energy 
transfer between waves ; and the second is the generation of organized flow 
with vortical pattern with strong stability and coherence of form. A first ap- 
proximation is to break up the field (the electric potential (f)) in two distinct 
elements : vortices and random waves. We will suppose that there is a finite 
number N of vortices in a two-dimensional plasma and that these vortices 
have random position. The vortices are individually affected by the turbulent 
background. The turbulent background, in turn, is affected by the presence 
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of the structures at random positions. In addition the turbulent background 
has statistical properties generated by the nonlinear nature of the waves in- 
teractions, even at amplitudes below those necessary to condense vortices. 
Random growth and decay of modes at marginal stability is included as a 
drive with Gaussian statistics. 

Consider the field 0^°''**'=^'^ = 0^ of the discrete set of N vortices individ- 
ually represented by the electric potential 0^"-' {x, y) 

N 

ct>y{x,y) = Y,^f\x,y) (1) 

a=l 

interacting with a turbulent wave field 0""^^^* = 0. The total field is 

N 

^{x,y) = Y,€H^.y) + ^{^.y) (2) 

a=l 

and we want to determine the statistical properties of the field ip (x, y). We 
will construct an action functional and we will calculate the generating func- 
tional of the irreducible correlations (cumulants) of the fluctuating field (p. 
The action functional is expressed in terms of the field which we see as 
composed from vortices and turbulence, Eq.Q. For example, the two-point 
correlation is composed of four terms 

(^(pLp'j ortices ortices ^ _j_ ^^Vortices ^waves^ 

I ^^waves ^Vortices^ _|_ ^^waves ^wavess^ 

Our procedure consists in absorbing the two intermediate terms into the first 
and the last terms. This is done by calculating the auto-correlations of each 
component taking into account the presence of the other. In this operation 
the second and the third terms, although are identical in Eq.Q are regarded 
differently: the second term is seen as the contribution of the turbulent 
background to the auto-correlation of a gas of vortices and the third term 
is seen as the contribution of a set of vortices to the auto-correlations of a 
turbulent field. 

Therefore we consider that the action functional is composed of two dis- 
tinct parts: Sytp = the action for the system of vortices interacting with the 
random field; and S^v = the action for the random field interacting with the 
vortices. Then the statistical ensemble of realizations of the fluctuating field 
ip {x, y) will actually consists of the Cartesian product of two distinct parts. 
The generating functional will be 

Z = Zv<fiZ^v (4) 
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where each factor is calculated using the action defined above. Although 
the generating functional is factorized, which corresponds to splitting the 
statistical ensemble into two parts, these two parts are not independent since 
each will be calculated such as to include the effect of the field from the other 
subsystem. 

The interaction with an external current must be included in each of the 
action in order to calculate the correlations using functional derivatives. In 
the full action the current is introduced by adding a term in the Lagrangian 
density, J if, where (p is the total field, vortices plus random waves. When 
we separate the two components we have ip = 0^°''**'=^* -|- ^^a-^^s _ 
j^Vorticesj^j^waves _ rjj_^jg exprcssiou must bc inserted in the action functional 
for the full system, before the intermediate terms in Eq.(j2I) are absorbed into 
the first and the last. This means that the same current J will appear in the 
final expressions of the two generating functionals Eq.(j3} 

Z [J] = Zy^ [J] Z^y [J] (5) 

It will be shown below that the value of the field at a particular point can be 
obtained by functional derivation to the current at that point p = jj. We 
obtain the average of the fiuctuating field as 



J=0 



Z^y [J = 0] 6 J 



Zipv [J] 



(6) 

j=o 



The two-point correlation is obtained from the generating functional by ap- 
plying two times this functional derivative, with the current in two distinct 
points and taking finally the current to zero. 

We now explain how this will be done effectively. 

1.2 Outline of the present analytical approach 

1.2.1 Gas of vortices in the turbulent background 

In the case of the vortices interacting with random waves, the basic reasoning 
takes into consideration two distinct elements. 

The first is concerned with the statistical properties of a collection of 
discrete vortices with zero, or very short, range of interaction. At this stage 
of the problem the particular shape of the potential distribution in a vortex 
is not essential and will be simplified in those situations where only the 
positions of the centers are important. The statistical ensemble consists of 
the configurations of randomly placed vortices, seen as a dilute gas. 
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The second aspect is concerned with the interaction between one (generic) 
vortex and the random waves in the surrounding turbulence and calculates 
the effect on the form of the potential of the vortex. Instead of the exact 
solution we will have now a form resulting from the scattering due to the 
random perturbation produced by the turbulent background. 

The statistical properties result from the combination of the two elements, 
which will now be explained briefly. 

The first problem is very similar to the Coulomb gas in two dimensions 
or to any problem related to dilute gas of interacting particles. The partition 
function is calculated using the sum of the energy of the individual vortices 
and the energy of interaction. The first part can be calculated from the 
exact analytical solution of the differential equation of the model (Eq. ()10|) 
below). The second part contains the sum of the energies of pair interaction, 
a problem that is principially complicated in our case, i.e. in turbulent 
plasma. Our real vortices are neutral (they are simply a deformation, stable 
and regular, of the scalar potential of the velocity) and there should be no 
interaction of Coulomb or Ampere type. However there is an elastic medium 
between the vortices and any motion of one of them generates sound waves 
that may couple with the potential of other vortices, influencing their motion. 
This may suggest that the interaction is of the same nature as the Casimir 
effect but this is beyond our primary interest here. What we minimally need 
to include in our partition function is the effect of scattering of the vortices 
at close encounter, since we do not take into account in the present treatment 
the variation of the number of vortices due to merging (which would imply a 
chemical potential). There is an intrinsic spatial scale in all aspects related 
with drift-wave vortices and this is the sonic Larmor radius, altered by the 
combined effect of the diamagnetic and a uniform flow. We show below that 
the pair interaction contains a kernel with fast spatial decay, Kq {r/pg) (the 
modified Bessel function). There are several reasons to assume that this is a 
physically correct choice (see jH])- 

The second aspect of the problem of vortices acted upon by turbulence is 
related to the direct modification of the exact vortex by the random waves 
around it. The case of a single vortex interacting with random waves has 
been treated in the papers [2], The starting point is the exact vortex- 
type solution of the differential equation we investigate. This is of course a 
non-random object and when it is isolated the statistical ensemble is trivial. 
However its interaction with random waves of the background turbulence 
makes it also a fluctuating object. We calculate the statistical properties 
of the ensemble of states of the fluctuating potential corresponding to the 
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shapes of the vortex-type solution. 

The action functional is extremum at the vortex solution, whose explicit 
expression we will use. A purely turbulent field also realizes the extremum of 
the action, but our separation between structures and random waves amounts 
to an approximative representation of the turbulent field, as a complimentary 
part to the structures. All configurations consisting of structures with ran- 
domly deformed, fluctuating shapes, together with random waves are a very 
good approximation of reality and must be found in close proximity of the 
extremum of the action. Then we have to explore the functional space of the 
system's configuration in the neighborhood of the structure (the extremum) 
trying to include as many as possible nearby configurations in the calculation 
of the generating functional. This will include into the generating functional 
the turbulent field, besides the structure. In few words, the idea is that the 
structure and the random waves, although very different in geometry, share 
a common property: they obtain or are very close to the extremum of the 
action functional. The technical procedure will consists of expanding the 
action functional to second order around the structure and to integrate over 
the space of configurations. This will automatically exclude bad approxima- 
tions, i.e. the states which are too far from realizing the extremum of the 
action since for them the Boltzmann weight is exponentially small. 

We conclude by noting that this is the standard semiclassical treatment 
0. In this sense it is similar to the treatment of vortex statistics of the 
Abelian-Higgs model of superfluids [S] and of many other systems. 

1.2.2 The turbulent field influenced by random vortices 

For the turbulent field interacting with vortices, our approach combines two 
distinct elements. The first is the inclusion of the vortices as a random per- 
turbation in the generating functional of the turbulent waves. The second is 
a perturbative treatment of the intrinsic nonlinearity of the turbulent waves, 
already modified by the inclusion of the effect of the random vortices. 

In the first part we follow the similar approaches as for the electron con- 
duction in the presence of random impurities, or as for flexible polymers 
in porous media. For the case where the vortices have uncorrelated ran- 
dom positions and take at random positive or negative amphtudes of equal 
magnitude we find that the problem is mapped onto the sine-Gordon model 
(actually sink). For our model equation this will amount to a renormal- 
ization of the coefficient which plays the role of a physical "mass" of the 
turbulent field ^01, |HH|, or, in other terms to a shift of the spatial scale from 
Pg (1 — Vd/u)~^^'^ to higher values. 

The second part consists of a systematic perturbative treatment of the 
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nonlinearity in order to get, as much as possible, a correct representation 
of the nonhnear content of the field of the turbulent waves (this means the 
nonlinear interaction and nonlinear energy transfer between low amplitude 
random waves [Hj , [?] , jl2I , llHj ) • The nonlinearity is included in a functional 
perturbative treatment where the turbulent plasma is driven by random rise 
and decay of modes at marginal stability. This induces a diffusive behavior 
of the turbulent field, at lowest order. We develop the treatment to one-loop, 
which means of order two in the strength of the nonlinearity. 

One may inquire if this is not in contradiction with the separation op- 
erated at the beginning, where the most characteristic aspect of the nonlin- 
earity, the generation of structure, is treated separately. However, it is clear 
that there is no danger of overlapping and double counting of the nonlinear 
effects. Since a statistical perturbative treatment is essentially an expansion 
in a parameter representing the departure from Gaussian statistics we can 
only hope to include higher cumulants beyond the second one (which means 
Gaussianity) but only few of them are accessible to effective calculation. Or, 
any structure needs a very large number of cumulants since, by definition, is 
an almost coherent field. It is illusory to try to capture the structure using a 
perturbative treatment starting with a state of equilibrium {i.e. no perturba- 
tion or, alternatively, a Gaussian collection of linear waves). Since any term 
of the perturbation can be represented by a Feynman diagram, we face the 
well known problem that the proliferation of diagrams at high orders leads 
to an effectively intractable problem. This is actually one reason for the use 
of the semiclassical methods. 

2 The physical model 
2.1 The equation 

The model of the ion drift instability in magnetically confined plasmas can 
be formulated using the fluid equations of continuity and momentum con- 
servation for electrons and for ions. It has been shown by a multiple space 
time scale analysis that the dynamics is dominated by two nonlinearities: 
the Charney-Hasegawa-Mima type, or vectorial nonlinearity, generated by 
the ion polarization drift and the Korteweg-De Vries, or scalar nonlinearity, 
related to the space variation of the density gradient length. The former 
is of high differential degree and is dominant at small spatial scales, of the 
order of few sonic Larmor radii p^. The latter is dominant at "mesoscopic" 
spatial scales, of the order of ^J'p^L^. The numerical studies [?], the fluid- 
tank experiments [?] and the multiple space-time scale analytical analysis 
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show that the scalar nonhnearity becomes dominant at late regimes in 
the statistical stationarity of the drift wave turbulence. However the possible 
manifestation of the two types of nonhnearity rises the problem of structural 
stability of either regime where only one of these nonhnearity is considered 
dominant: inclusion of the other strongly changes the behavior of the system. 
It is then pertinent to consider that the turbulence generated in a realistic 
regime may include manifestations of both types. The turbulence is domi- 
nated by the larger scales sustained by the scalar nonhnearity (described by 
the Flier 1-Petviashvili equation ^3]) together with robust vortices generated 
at the scales of few p^, typical the CHM {i.e. vectorial) nonhnearity [?]. 

When the scalar nonhnearity is prevailing [?] the equation has the 
form 

A Y72\9f d(f dip 

(l-Vj^ + ^^*^-^^.T^^ = (7) 
In a moving frame and restricting to stationarity we obtain 

Viv? - acp - /3ip'^ = (8) 
The physical parameters are [0], [Ej, [ini 

where = Cg/^i, Cg = {Te/mi)^^"^ and the potential is scaled as — > e^/Tg. 
Here L„ and Lt are respectively the gradient lengths of the density and 
temperature. The velocity is the diamagnetic velocity Vd = PgCg/Ln- The 
condition for the validity of this equation are: {kxPs) [kp^f <^ rj^pjLn, 
where t]^ = Ln/ Lt^. The coefficients a and (5 have the dimension {length)~'^ . 
This form will be used below. 



2.2 The structures 

The exact solution of the equation is 



Vd 



X sec h^ 



2^ {^-'iT' (y-yo-^t) 



(10) 



where the velocity is restricted to the intervals u > Vd or m < 0. In Ref . |14j 
the radial extension of the solution is estimated as: (Ax)^ ~ Ps-^n- In our 
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work we shall assume that u is close to Vd , u > Vd (i.e. the structures have 
small amplitudes). The monopolar vortex in this regime is discussed in Ref. 
[?]. For asymptotic form of the CHM equation see Ref.jH]. We will adopt 
the one-dimensional section of the solution {i.e. Eq. ljlOj) ) when we calculate 
the eigenmodes of the determinant of the second functional derivative of the 
action. In conclusion, we consider coherent structures which are monopolar 
vortices of both signs of vorticity, with equal magnitudes and with random 
positions in plane. 

3 Statistical analysis of the physical system 
3.1 General functional framework 

The general method for constructing the action functional for a classical 
stochastic system is described by Martin, Siggia and Rose (MSR) [22], and 
in path integral formalism, by Jensen |2H1- Two reviews by Krommes are 
very useful references on this point jSHj, [HO]- The functional method has 
been applied in several concrete problems and references may be consulted 
for details |S2], IMI; [2|- We here review few elementary procedures 
(see ESI). 

Consider a differential equation F [(f)] =0 whose solution is 0^. The un- 
known function belongs to a space of functions (f){x,y). We want to select 
from this space of functions precisely the one that is the solution of the dif- 
ferential equation and for this we can use the functional Dirac 6 : 5{(j) — (j)^). 
This can be represented as a product of usual 5 functions in every point of 
space 

N 

K4> ' ct>^) = X{nct>{xu) - ct>^ {xu)] (11) 

k=l 

Any operation that will be done on a functional of </> can be now particularized 
to the solution (f)^ by simply inserting this Dirac functional. For example the 
calculation of a functional VL (0) at the function cj)^ can be done by a functional 
integration over the space of all functions, with insertion of this 5 functional. 
Using the Fourier representation of the ordinary Dirac functions and going 
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to the continuous limit we note the appearence of the dual function % 

6F 



5 [F (0)] 



6F 



X exp <i dxx (x) F [(f) (x)] 



121 



We will define the partition function as usual, by the functional integral 
of Boltzmann weights calculated on the base of the MSR action 

Z[J]^ Jv [0] V [x] exp I j dxdy {xF [0] + J^<f> + J^x) } (13) 

The functional integration takes into account the fluctuations of the phys- 
ical field and of its dual x- The "free-energy" functional is defined by 
exp {W [J]} — Z [J] from which the irreducible correlations (cumulants) are 
calculated by functional derivatives to J. 

(0 (x, y)) = exp {-W [J]} f exp {W [J]} (14) 

and similar for higher cumulants. The two-point irreducible correlation for 
the field is 



(0 {x,y)(j){x',y')) 

J V[<f>]V[x]<f>{x,y)<f>{x',y') 



X 



exp I J dxdy {xF [0] + J^4> + J^x) 



j=o 



(15) 



Z[J ^0] 5J^(x,y)5J^(x',y')' 

This is the general analytical instrument that will be used in the following 
calculations. The calculations from this work are presented in greater detail 
in Ref. [?]. 



4 Vortices with random positions 

4.1 The discrete set of vortices in plane 

The action of the discrete set of vortices is determined by the sum of the 
actions of the individual vortices plus a part that results from the interaction 
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between them . The first part is simply the time integration of the energy, 
i.e. (since time factorizes) the space integration of the expression of the 
product of the static potentials associated with a single vortex, Eq. (fTO|l . and 
to its dual, which in the end simply means the square of the wave-form of 
the vortex potential. This quantity is repeated for each of the vortices. 
The partition function is 



with the following meaning. 

The first factor simply takes into account independent vortices with 
arbitrary positions in plane and expresses the fact that this part of the parti- 
tion function results from a Cartesian product of the statistical ensembles, 
one for each vortex. The factor A^! takes into account the permutation sym- 
metry. The generating functional for a static vortex with structure given by 
the interaction with random waves, is calculated in Eq. ()65|) below. We have 
Zy^ = Zy^ [J] where we have indicated by the index V(p that the partition 
function of the vortex includes the interaction vortex-turbulence, and that 
the expression depends on the external current J, and will contribute to any 
correlation that we will obtain by functional derivations at J. 

The sum is over the set of configurations {a} characterized by random 
choices of positive and negative vortices. 

The integrations over the positions of the centers R*-"^ of the vortices, 
a = 1,N, express the fact that we allow arbitrary positions in plane, with 
equal probability. Each integral is normalized with the area of the physically 
interesting two-dimensional region of the plane, A. 

The exponent of the Boltzmann weight contains action resulting from 
the interaction between vortices. We expect that for a dilute gas of vortices, 
where the distance between the centers R*-''^ — R^**^ is much larger than the 
core diameter d, R^"^ — R*^''^ 3> c? , a,b = 1, N , the interaction is very weak. 
In order to describe the interaction between vortices we start from the well- 
known alternative model of the drift waves sustained by the ion polarization 
drift nonlinearity [TSj 12^1 • In this model it is considered a set of A^^ 
point-like vortices of strength Ui interacting in plane by a short range poten- 
tial expressed as the function Kq (modified Bessel function) of the relative 




(16) 



12 



distance between vortices. The potential (f)^ in a point R is a sum of con- 
tributions from all the A^^ point- vortices (jf (R) = Xli^i ^j-^o (P7^ l^- ~ R-jI) 
and the equations of motion dUi/dt = —'V(ff x (where is the versor 
perpendicular on the plane). The distribution of vorticity of the physical 
system (in particular the quasi- coherent vortical structures) represents spa- 
tial variations of density of these point-like vortices. The interaction between 
the physical vortices will result from the interaction between the point-like 
vortices, taking into account the density of these objects. The energy of 
interaction is 



H = 2^ 2^00,00 J Ko {p;' |R. - RjI) 



i>j 

called the Kirchhoff function. The range of spatial decay of the interaction 
is the Larmor sonic radius p^, which however may be modified to an effective 
Larmor radius, in the presence of gradients and flow. When we approach the 
continuum limit N^^ — > oo the envelope of the density becomes the physical 
vorticity uj {x, y) which, for this stage of the problem is sufficient to be con- 
sidered as highly concentrated in the cores of the physical vortices and almost 
vanishing in the rest. Taking the elementary point-vortices of equal strength 
ujj = cuq we have that each physical vortex is an integer multiple A^*^"-* of 
this quantity. Now we will associate with each physical vortex a continuous 
function, i.e. its vorticity defined on the whole plane, p!^^ {x,y), which is, as 
said, concentrated in {x, y) 

p(r)(a:,2/) = iV(«W(R-R^'^^) (17) 

Then the energy is [S] 

TV N 



H = / dxdy / dx'dy' 

a=l b=l 



a>b 



I "D (a) -p (b) I 

xp(;) (x, y) K, I I 1 1 (x', y') (1^ 



Due to Eq. ()17|) the interaction energy is only the interaction between the 
centers R^""^ and R^^^ of the vortices. The summation proceeds by grouping 
the point-vortices into physical vortices, then assuming that these (for only 
this stage of the problem) have 5-function shape and finally formally replacing 
this 5 with a continuous distribution pL"'' [x, y). In this operation a number of 
infinities arise from the energy of the interaction of the point-vortices which 
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are grouped into one physical vortex, since the relative distances are zero for 
them. This singular part can be removed since it does not participate to the 
functional variations induced by the "external excitation" current J. 

For simplification of the computation we now only consider physical vor- 
tices of equal amplitude (positive or negative) and then the vorticity distri- 
bution pif^ (x, y) has amplitude a;„ = pujQ (p is an integer), multiplied by the 
integer r?/") which can take the values ±1 for positive or negative vorticity. 
WehavciV^") = pn^'^l The sum over the physical vortices' positions suggests 
to define a formal unique function of vorticity {x, y) for all the N physical 
vortices 

N N 
a=l a=l 

Further, the energy is normalized with a constant dimensional factor. The 
interaction part can be rewritten 



exp 



N N 



-^flY^Y. ^^^^ / 



a>b 



exp 



r, 4 2 / ^^^y [ dx'dy'p^{x,y)G (R-R') p^{x',y') 



(19) 



where G is the kernel of interaction, 

G(R-R') ^ i-5;5;if„(p.|R(")-R(')|) 

xd (^r-r(")) s (r'-r^''^^ 



(20) 



The differential equation for Ko is (A - l/p^) Ko (r/p^) = -27:5 (r) . This 
helps to replace the interaction part with a Gaussian functional integral, by 
introducing an auxiliary field ip 



exp 



■-^ / dxdy f dx'dy'p,{R)G{R-R')p^{R') 



= Pi^ j ^^Mcxpj-^ j dxdy 

dxdyp^ (R) %Ij (R) 



pI 



X exp 



pIi^v 



(21) 
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with pi a normalization constant. We make a change of variable in the 
functional integration 2ttiI) —>■ x (this also changes the normalization constant 
Pi p) and return to the partition function Eq. (fT6|) 



p- 



V [x] exp 



dxdy 



yxf + 



X — r I ^1/ 



(0) 



{a} \a=\ 



exp 



:22i 



In the last factor we note that in the sum each term consists of two contribu- 
tions, corresponding to positive and negative vorticity, n^"'^ = ±1, and they 
are weighted with the same factor, 1/2 



{q} a=l 



7(0) 



N 



' ^ r 1 1 

n / jdR- (exp [ix (R)] + exp [-tx (R)]) 

a=l 



(23) 



In the last line we have removed the upper index (a) since all factors in 
the product are now identical. For a fixed number of vortices the parti- 
tion function Zy^ (with nonvanishing contribution to the derivatives to J) 
is decoupled from the other factors and will provide iV-times the same con- 
tribution. The other factors, i.e. the functional integral that contains the 
interaction between the vortices can only appear in the final answer as a con- 
stant, multiplying contributions coming from Zy^ [J]. When is arbitrary 
the partition function must also include a sum over terms each corresponding 
to a number of vortices 



Af=0 



exp 



.(0) 



N 



7(0) 
A 



jdK^ (exp [ix (R)] + exp [-tx (R)]) 



dxdy cos [x {x, y)] 



N 



(24) 



Then the partition function becomes 

1 



Zy=p-'Jv[x]e^pL^^J 



dxdy 



A 



Zy^ COS [x {x,y)] 
(25) 
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In this expression the quantity Zy is a functional integral over the space of 
fields (x, y) representing a single vortex. In the absence of the background 
turbulence the field is a deterministic quantity [Eq.lfim)] and the 

statistical ensemble is trivially composed of one element. The interaction 
with the background turbulence induce a fluctuating form and the statistical 
properties can be obtained from Zy^ . In other words, Zy includes two sources 
of fluctuations: one is the fluctuation of the field of vorticity due to the 
random positions in plane of the vortices (the gas of vortices) and the other 
is the fluctuation of the shape of the generic vortex due to interaction with 
the turbulent background. 

If we want to use the expression as a generating functional for correla- 
tion we must be able to drop into the functional integral the field representing 
the vortices, i.e. 

c^y (R) = 4^^:^ (x, y) (26) 

a 

as in Eq. ()15|) . An external excitation by the current J will produce a change 
in Zy from the change of the vorticity of the gas of vortices and from the 
change of Zy\ We have 



5Z,, 



6 J {x,y) 



+ 



Vort 



SZy 6Zf 

6ZpSJix,y) 



(27) 



{(py {x,y) (t)v {x',y')) 



Zy [j = 0] 
V 



5^Z 



V 



6J{x,y)6J {x',y') 



:28) 



vort 



6'Z, 



5Z. 



(0) 



5Z. 



(0) 



m\'^5J {x,y)6J {x',y') 



6 Z 



V 



6Z 



V 



6'Z^^^ 



SZfSJ{x,y)6J{x',y') 



The formulas are taken at J = 0. The first terms in these equations are 
related with the fiuctuations of the vorticity of the gas of vortices as a con- 
tinuous version of the discrete set of physical vortices with arbitrary positions 
in plane. The other terms are related to the fiuctuation of the shape of a 
vortex and in order to calculate these contribution we need the explicit ex- 
pression of Zy, as a functional of Zy\ Further, we will need the detailed 
expression of Zy^ [J] and this will be calculated in the next Section. 

In order to obtain the contribution from the fiuctuation of the gas of 
vortices (the first terms in Eqs. ()27|1 and (j^H|) ). we introduce a new term in 



16 



the action, consisting of the interaction between the vorticity distribution 
(R) and an external current, 

i— I dxdy[p^{x,y) J^{x,y)] (29) 

(the factor i is introduced for compatibihty with Eq. ()21|) ). This current 
is an external excitation applied on the field of the vorticity and not on the 
field of potential we would need in the Eq.(jni). We may assume that 

there is a connection between the current J^^ and the current J (which acts 
on the field (py) but there is no need to specify this relation. Indeed, the 
Eq.® shows that at the end both currents should be taken zero. 
The last line of Eq. ()21|) transforms as follows 



exp 



i— I dxdyp^(R)ij (R) 



exp \i— [ dxdyp^ (R) [^/j (R) + (R)] 

I i^v J 



(30) 

All the calculations following Eq. 1)211) are repeated without modification but 
in the last term, instead of the function x {^^ u) "we will have 

cos [x (x, y)] cos [x (x, y) + (x, y)] (31) 

since this was the term which resulted from the presence of the function 
(R) in the Eq.(j2^. Making the change of variable in the functional inte- 
gration (of Jacobian 1) 

X ^ X - ^ (32) 
the integrand in the action is expressed as 

{ [V (x - J.)f + j,{x-J.f- ^4°^ COS [x (x, y)] I (33) 



The two-point correlation of the field of the vorticity fluctuations can be 
calculated from 



ujI ' ' Zv[Juj = 0]i6J^{x,y)i6J^{x',y') 

The squares in the action Eq. fjHHj) are expanded 



(34) 



^ ''[Vx)' + ^x'-^4°^cos[x(x,y)] (35) 



87r2 p2- ^ 

+ (V JJ' + ^ 

Ps 

-2 (Vx) (VJJ - ^xJ. 

Ps 
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In the part of the action that depends on we make integrations by parts 



(V JJ' + -,Jl-2 (Vx) (V JJ - -xJ. 



Pt 



(36) 



-J^ (AJ^) + \jl + 2^ (Ax) - 



The first hne at the exponent in Eq. ()35|l does not contain the current and in 
the following functional derivations to we will temporary omit it. Consider 
the apphcation of the first operator of derivation to i (a;, y) 



i6J^ {x,y) 



exp 



87r2 



dxdy 



-J^ (AJ^) + ^Jl 



1 

7s 



+2J^ (Ax) - ^X^ 

Ps 



(37) 



exp{...} 



xp: 



-1^ 



-2AJ^ + 4^ + 2Ax-4x 

Ps Ps / 



Every derivation to the current suppresses a space integration and in 
consequence the result is multiplied with factors which render the space 
integral dimensionless. The subscript shows that the functions inside the 
bracket are calculated in the point {x,y). 

The second operator of derivation is now applied on Eq. ()37j) 



6 

■cj I , A exp 



1 



87r2 



dxdy 



-J^ (AJ^) + — J2 



+2J^ (Ax) - ^X^ 

Pi 



1 2(-i. 



^7^^ 8vr2 

1 .hi) 
87r2 



I 2(-l) 

X exp {...} 

II 2(-l) 



+ --Ps 



I I- ' 87r2 



X exp {...} 



-2AJ^ + 4^ + 2Ax-4x 

Ps Ps J {x,y) 

-2AJ^ + 4^ + 2Ax - 4x , 

Ps Ps / 

-2AJ^ + 4^ + 2Ax-4x 

Ps Ps / 



-2A5 {x -x',y- y') + —5 {x - x\ y - y) 

Ps 



{31 
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At this moment we can take J = 0. The exponentials in the Eq ()38|) are equal 
to 1. The auto-correlation is 



out 



(39) 



:Vxf + 4x^-^4°^ cos [X. 

Ps ^ 



-2A6 {x -x',y- y') + —5 {x - x', y - y') 

Ps 



.2 \ 2 



Ps J {x,y) V Ps J {x',y') 



The last two lines in Eq. ljH^ (the curly bracket) arise from derivation. The 
normalization gives by definition 



Zv[J^ = 0] = (1) 



Zy COS [x (x, y)] 



The functional integration does not affect the first two terms in the curly 
bracket in Eq. ()39|l . The functional integral of the last term represents the 
average of the product of two functions x- 



P. 



{Pu ix,y)Pu W,y')) 



47r2 



-4 



-A5 (x -x',y - y') + —5{x- x\ y - y') 

Ps 



P. 



2 \ 2 



A 



1 

Ps / {x,y) 



A 



1 

PsJ (x',y') 



(x {^,y)x{^' ,y')) 



It is now easier if we use the Fourier representation of the fields, which 
we denote by the symbol~ Writting x = {x,y), 

\{pu>{x,y) Pu>{x\y')) = \ I c?kexp(ik-x) / c?k'exp(ik'-x') (p<^(k)p^(k')) 

J J 
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We have 



^ j dkexp (ik ■ x) j cik'exp (ik' • x) (p^ (k) (k')) 
y rfkexp [ik- (x - x')] (^P + -^^ 



2 \ 2 



/.kexp(.k.x)/.k'exp(.k'.x' 



x(A:2 + - A:'2 + -)(x(k)x(k')) 



47r2 

We can take a fixed reference point 

x' = a 
X = a + X — x' 

and write 



\ f dkexp [ik- (x - x')] / dk' exp [z (k' + k) ■ a] (p, (k) (k')^40) 
j dkexp [ik- (x - x')] [k"^ + ^ 



^2 



p: 

47]-2 



2 \ 2 
P. ^ 



4772 



y dkexp [ik- (x - x')] j rfk'exp [i (k' + k) - a] 



Tlie parameter a lias no particular role : non of our assumption has imposed 
a nonuniformity of the statistical properties on the plane. Therefore we can 
integrate Eq. ()40|) over the position a, i.e. on the plane 



1 



da. 



Obviously, this will produce in the left hand side a function 6 

6 (k' + k) 

after which the integration over the second wavenumber, k' , will impose 

k' = -k 



20 



For the first term in the right hand side, the integration over a will have no 
effect. For the second term the effect is the same as in the left hand side, i.e. 
we have k' = — k. we will now replace x — x' by x and obtain 

■\ I dkexp(ik-x) (p^(k)p^(-k)) 



47]-2 



dk exp (zk ■ x) { k"^ -\ — - 

pi 



jdk exp (.k . x) + ' (X (k) X (-k)) 



In physical space the correlation also reflects the statistical uniformity. 

If _ _ 1 

— / dkexp {ik ■ x) (p^ (k) p^ (-k)) = — (p^ (x) p^ (0)) 



out 



The equation is 



XiPu. (k) P^(-k)) 



(41) 



P 



V 

2 / 1 
2 ^ 



47^2 



pI 



47r2 



Pt 



+ (X(k)x(-k)) 



The second term in the bracket of Eq. ()4H) contains the two-point correlation 
of the function x and can be obtained by explicit calculation of the functional 
integration in Eq.(??). The same analytical problem as for the explicit cal- 
culation of Zy [J = 0] and (xx)k = ix (k) X (~k)) will appear later (for the 
turbulence scattered by the random vortices) and there we will give some 
details of calculation. At this moment few explanations are sufficient. For 
small amplitude of the auxiliary field x the function cos is approximated with 
its first two terms ^ 

cosx^l-y (42) 

The constant 1 is only a shift of the action. However it leads to a term that 
depends on Zy\ which is integrated in the exponential over all volume i.e. 
the area A on the plane. In detail, replacing Eq. ()42j) in Eq.(??) 



Zv [J. = 0] 



-1 



P 



X exp 



X exp 



»7r^ 
1 

'8^ 



dxdy ^- 
dxdy 



(43) 



yxf + -x' + 

Pi 



Svr^ „(o)\ X'' 



A 
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The first factor can be taken outside tlie functional integration 



exp 



exp 



dxdy I - 



(0) 



(44) 



7(0) 



Since it is determined by the non-interacting vortices, it must exist even if we 
would neglect completely the interaction between the vortices taking x ^ 0. 
The rest of the Eg. ()43j) is the Gaussian functional integral 



P 



-1 



(45) 



X exp 
X exp 



87r2 
1 

'8^ 



dxdy 



yxf + -a' + 



yio)\ X' 



A 



(46) 



dxdyx {x, y) 



-A 



An' 



(0) 



x{x,y) 



where 



1 47r 



7(0) 



The determinant can be calculated explicitely, by the product of the eigen- 
values of the operator. This product (besides an infinite factor that will 
disappear) is convergent. However we keep this formal expression 



Zv = Zv [J^ = 0] = q exp [det (- A + a^)] 



,2M-l/2 



(47) 



The second term in the bracket of Eq. ()41|) contains the two-point corre- 
lation of the function x and can be obtained by explicit calculation of the 
functional integration in Eq. ()39|) . The same analytical problem as for the 
explicit calculation of Zy [Ju = 0] (Eq.(??)) and (xx)k = ix (k) X (~k)) will 
appear later (for the turbulence scattered by the random vortices) and there 
we will give some details of calculation. At this moment few explanations 
are sufficient. For small amplitude of the auxiliary field x the function cos is 
approximated with its first two terms 



cos X ~ 1 — 



X^ 



(48) 



22 



The constant 1 is only a shift of the action. However it leads to a term that 
depends on Zy\ which is integrated in the exponential over all volume i.e. 
the area A on the plane. In detail, replacing Eq. ()42|l in Eq.(??) 



Zy [J. = 0] 



1 



1 



dxdy 



X exp 



dxdy 



(0)" 



The first factor can be taken outside the functional integration and is exp 

Since it is determined by the non-interacting vortices, it must exist even if we 
would neglect completely the interaction between the vortices taking x — 0. 
The rest of the Eq.ipHj) is the Gaussian functional integral. Introducing the 
notation = l/pl + An'^Z^^ /A 



Zv = 0] = gexp zf [det (-A + a")] 



-1/2 



(49) 



and g is a constant. The determinant can be calculated explicitely, by the 
product of the eigenvalues of the operator. We need this explicit expression 
because we need the functional dependence of Zy = Zy [J^^ = 0] on Zy^ 
as results from Eq. ()28|) . As will become clear later, the factor with the 
determinant, which comes from the influence of the fluctuating shape of a 
vortex on the correlations of the vorticity of a gas of vortices in the plane, is 
affected by a factor p^/A, which is small compared with the exponential in 
Eq. ()47|) . Therefore we calculate in a one dimensional cartezian approximation 
the eigenvalues, instead of a cylindrical problem. We have to solve 



d^ 



dx"^ 



on an interval L. The eigenvalues are = {2TTn/Ly + ^ where n IS an 
integer, and we obtain 



det(-A + a2) = Y[xi = Y[[{2nn/Lf + a^] 



n=l 

oo 



n=l 



1 + 



n {27rn/Lf J] 

n=l n=l 

sinh (aL/2) , 

' n=l 



a^Ly{27ry 



(50) 
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The infinite product is eliminated since we always use the ratios of Zy [J^ 
and Zy [ = 0] . We also take L = a/A and obtain 



sinh 



Zy = q exp 



.(0) 



2^(0) 
V 



1/2 



/2 



-1/2 



(Ap72 + 4vr24°))'^'/2 



(51) 



The quantity Ap^ ^ is very large and an approximation is possible 



sinh 



[Apj' + Att'Z 



2 7(0) 
V 



1/2 



/2 



(^Ap-^ + An^Z^y' /2 



exp 



exp (p;'VA/2) / 

— ^ 7= ^exp 

p-^VA \VAp 



Ps'VA 



-1 



2 



-2 



Api 



(52) 



The first factor is large but constant and can absorbed into the coefficient q. 
We get in Eq.(jSH) 



gexp 



gexp 



7(0) 



exp 



Zf I 1 



2v^p-i 



-1/2 



-2 



2v^p 



-1 



1 + 



Ap 



-2 
s 



We can neglect the second term in the first exponential 



Zy = [J^ = 0] = g 1 



Aps 



exp 



7(0)' 



(53) 



The second term in Eq. ()4H) . i.e. the auto-correlation of x k-space, 
may be calculated starting from the real-space correlation 



Zy [J = 0] 



X exp 



p ^ exp 



dxdyx (x, y) (- A a^) x {x, y) 
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As usual we return to the form of Eq. ()45|) . and only for this step, we insert 
an external current Jg {x, y) interacting with x- The auxilliary functional is 
denoted [Jg] 



Zv [J = 0] 
X exp 
1 



p exp 



AO) 



dxdy [x (x, y) (-A + a^) x (x, y) + Jex] 



Zv [J = 0] 



p ^ exp 



X exp 



dxdycj) (a;, y) (—A + a^) (x, y) 



j dxdy^Je{x,y){-A + a^) ^ Je{x,y)^ 

To obtain the above equation we have made a change of variables X = 
X + i (—A + a^) ^ Je of Jacobian 1. The functional integration over can 
now be carried out and the rest of the expression at the exponent appears in 
a factor 

~ [det (-A + a^)] exp i^--^ j dxdy^Je{x,y) {-A + a^) ^ Je{x,y) 

where the symbol ~ means that there also result constant factors. But these 
are the same as those contained in the factor q introduced in the Eq. ljTfj) . 
We then have 



Z^ [Je. 



1 



p exp 



Zv [J = 0] 

X exp I - ^ / ^^^^ [x{x,y) {-A + a'^) x{x,y) + J^x] 
p^^expj-^ j dxdy^Je{x,y) {-A + a^)'^ Je {x,y)^ 
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where we have taken into account Eq. ()47p . The correlation is 

1 S^Z^ [Je] 



Zx [Je = 0] SJe {x, y) 6Je {x' , y') 



Je=0 



1 



6 



Stt 



^2{-A + a') 'Je{x,y)]exp{...} 



J Je = 



[^e = 0] 



P 



-1 



47r^ 



(-A + a^) '5(x-x')exp{...} 



;2{-A + a') 'Je{x,y) 



2(-A + a2) ' J,{x',y')]exp{...} 



in' 



.(0) 



(the factor exp 
have 

{x{x,y)x{x',y')) 

or 



has not been written since it disappears). Finally we 



47r2 



47r2' 



47r2 



.(0) 



(54) 



(55) 



Introducing this in Eq. ()41|) the correlation of the field of the scalar potential 
will have from this a contribution (with normalization) 



I I \ vort 

V0y)k 



1 + 



(56) 



X 



47^2 



2 ^2 



1 + k'p^ 



1 + k^pl + 4°U7r2p2/^ 



which represents, as said, the first term of the Eg. ()28j) . 



The other terms in the correlation 



)(Eq. 



are due to the fluctu- 



ation of the form of the generic vortex from interaction with the background 
turbulence. We calculate, using Eq.(??), the derivatives 



5Z 



V 



and 



5^Z,, 



g 1 + 



2 7(0)~ 



+ 



V 



exp 



g 1 + 



27r2 vr^zf" 



Aps 



exp 



Now we have to calculate the derivatives of ^ to J as shown by the 
last two terms Eq.(j2Hl)- This part will be added after Zy^ is calculated, in 
the next subsection. 



7(0)' 



(57) 



(5^ 



J Je = 
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4.2 A single vortex interacting with a turbulent envi- 
ronment 



The partition function for a single vortex in interaction with turbulence is 
defined as 

=^'' j D [x] D [0] exp {Sv [x, 0]} (59) 

and the equation is the stationary form of the equation used in Ref. 0, 

F [(f)] = Wlif -aip- pif"^ (60) 

The density of Lagrangean 

^ [X {x, y) , (x, y)] = X (Viy? - a^p - (3^'^) (61) 

is obtained from the Martin-Siggia-Rose method for classical stochastic sys- 
tems. Then the action is 

Sv [x, = y dxdy^ [X {x, y) , (p (x, y)] (62) 

As usual we introduce the interaction with the external current J = (J^, J^). 
However since there is no need to calculate functional derivatives to Xj "we 
can only keep J = J^. For uniformity of notation in this paper, we will not 
use the factor i in front of the action in contrast with [2j. We have then 

4°) [J] =M-^ j D [x] D [0] exp I j dxdy [x {Vl^ - aip - ^i/) + J0] 

(63) 

We need the explicit expression of the functional integral Eq. (jU^ and this 
has been obtained in the references |2] and jS] • For convenience we will recall 
briefly the steps of the calculation restricting to the results we need in the 
present work. 

To calculate the generating functional of the vortex in the background 
turbulence we proceed in two steps: we solve the Euler-Lagrange equation 
for the action ()62|) obtaining the configuration of the system which extremises 
this action; further, we expand the action to second order in the fluctuations 
around this extremum (which will include the turbulent field) and integrate. 
The Euler-Lagrange equations have the solutions 

^jsi^^y) = fsi^^y) (64) 
Xjs{x,y) = -(Ps{x,y) + xj{x,y) 
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The first is the static form of the solution Eq. lfTIH) and does not depend on 
J. The dual function is — y?^ [x, y) plus a term resulting from the excitation 
by J in its equation. This additional term Xj (^) v) is calculated by the 
perturbation of the KdV soliton solution according to the modification of 
the Inverse Scattering Transform when an inhomogeneous term {i.e. J) is 
included. The action functional is calculated for these two functions 

Svs [J] = Sv [^s y)^-'P3 (x, y) + Xj (a;, y)] 

Then the first part of our calculation is Zy^ [J] ~ ex\i {Svs [J]]- Ex- 
panding the action around this extremum 



Sv [x, J] = Sv [Vjs, Xjs] + 



5^Sv [J] 



2 \ 5^5x 

we calculate the Gaussian integral and obtain 

5''Sv [J] 



[J] = AT-' exp {Svs [J]} 



det 



5(pSx 



-1/2 



If we can neglect the advection of vortices by large scale wave-like fluctu- 
ations, we can calculate the determinant since the product of the eigenvalues 
converges without the need for regularization. The result is 



4°) [J] = AT-i exp {Svs [J]} AB 



Where 



A = A[J] = 



B = B[.J] 



(3/2 



sinh [13/2) 
a/2 



1/4 



1/2 



(65) 

(66) 
(67) 



sin (a/2)_ 

The eigenvalue problem depends functionally on Xj i.^, y) which implies that 
/? and a depend on J. Their expressions can be found in 0. With those 
detailed formulas we can consider that we have the necessary knowledge to 
proceed to the calculation of the functional derivatives of Zy"* [J] at J, using 
Eq.dHHI). 



zP [J = 0] 



SZ!^^ [J] _ 6Svs [J] ISA 1 6B 
6J ~ 6J ^ aJj ^ bJj 



(68) 
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For the present problem it is sufficient to take the ffist term as Eq.()10|) 



.J - <Ps (69) 

j=o 



The next two terms in Eq. ()68j) represent the averaged, systematic modifica- 
tion of the shape of the field around the vortex due to the mutual interaction. 
They will be equally neglected, assuming that the main effect is contained in 
the dispersion of the fiuctuations of the shape of the vortex interacting with 
the random field (which actually is our main concern here). Then 

' ^^i-"'!^'... (70) 



Zf [J = 0] w 
The second derivative to the excitations in two points yi and y2 is 

1 [J] 



(71) 



ZP [J = 0]SJiy2)5Jiy,) 
6Svs[J]6Svs[J] S'Sy^J] 
SJ{y2) ^J{yi) SJ{y2)6J{yr) 

1 6 A 6Svs [J] 1 6B 6Svs [J] 
^A6J{y2) 6J{y,) + B6J{y2) 6J{y^) 

1 SA 6Svs [J] 1 6B 6Svs [J] 
'^A6J{y,) 6J{y2) ^ BSJ{yi) 6J{y2) 

1 SA 1 SB I SA 1 SB 
'^ASJ{y,)BSJ{y2) ^ ASJJ^BSJi^ 

1 S^A 1 S^B 

^ASJ{y2)SJ{y,) + B SJ {y2) SJ {y,) 

Since we have assumed as an acceptable approximation to neglect the av- 
eraged change produced by the turbulence on the soliton shape the second 
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term in the RHS is zero. For the first term we use Eq. ()69|) . We have 



1 



[J] 



Zf [J = 0] 5J(l/2)5J(i/i) 



0s (1/2) 0s (Z/l) + 

+0S (1/1 



+0S (1/; 



5J (2/2! 
5 



5 



(^2 

5 



■In 5 



5J (yi) 5J {yx) 

51nA (51n5 51nA 51n5 

(5J(yi)5J(y2) ^5J(l/2)5J(l/i) 

1 1 " 

+ 



In 5 



A5j(y2)5J(i/i) 5 5j(i/2)5J(i/i; 



The expressions are comphcated (see the Appendix of Ref.[3 ) and some 
numerical calculation of these expression is unavoidable. For small amplitude 
of the turbulent field the expression can be rewritten 

(0) 



[J = 0] 5^ (2/2)5^(2/1 
0s (2/2) 0. (2/1) [1 + / (y)] 



J=0 



i.e. in a form that expresses the fact that the two-point correlation is ba- 
sically the auto-correlation of the potential of the exact soliton modified by 
a function / which collects the contributions from the interaction with the 
random field. In fc-space we have 



S'zf [J] 



Zf \J = 0] (2/2) SJ {y. 



(k) 0, (-k) [1 + / (k)] 



J=0,k 



At the limit where we do not expand the action to include configurations 
resulting from the interaction vortex-turbulence, we have a — and /9 ^ 
and it results A = B = 1. In this case / = 0. 

The detailed expressions of these terms are given in the paper |2] , • For 
the purpose of comparisons we will express the spectrum as 



)k = 5(k) [l + /(k)] 



(72) 
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where 0q is amplitude of a vortex, S (k) has been derived by Meiss and 
Horton [H] 



and / (k) is function that is the correction to the Fourier transform of the 
squared secant-hyp erbohc, produced by the turbulent waves. 

Before proceeding further with the calculations based on the Eq. fl72|) we 
need to discuss the formal term / (k). Since this term represents the differ- 
ence from the simple isolated vortex to the vortex perturbed by turbulence, 
one would like to have a quantitative connection between the amplitude of 
this term and at least two elements characterizing the background turbulence: 
(1) the amplitude and (2) the spectrum. 

In the way we have conducted the calculations of the generating functional 
Eqs. (j65|) . (j66|l . (j67j) the new terms in the expression of the auto-correlation 
due to the factors A and B are expressed in real space, not in Fourier space. 
They are obtained from the product of eigenvalues of the operator repre- 
senting the second order functional derivative of the action, i.e. they are 
connected with the geometry of the function space around the exact, vortex, 
nonlinear solution. The determinant of the operator S'^Sy/ {S^^x) may be 
seen as a volume in the function space, centered on the vortex solution. The 
inverse of any eigenvalue gives an idea of the extension along a particular di- 
rection (eigenfunction) in function space. When an eigenvalue is very small, 
the operator almost vanishes on functions along that direction. At the limit 
this is a zero mode and corresponds to a translational symmetry of the phys- 
ical system along that direction. The correlations depend on the sensitivity 
of this volume (the product of the eigenvalues) on the excitation J applied on 
the system. The excitation is first manifested in the appearence of Xj {^^ v)- 
This one consists of a part that will modify the shape of the exact vortex 
plus the oscillating tail generated when a soliton is perturbed. The latter can 
be considered as a component of the background turbulence. The "propaga- 
tion" of the influence from an excitation J can be summarised symbolically 
in the chain : J — > Xj (x, y) eigenvalues of the operator S'^Sy/ {5(fSx) — > 
A and B (or cr and (3). The expressions of a and (3 can be found in jSj. 




(73) 



31 



Now we can return to the Eq.(j2HI)- Using Eq. ()70|) the second term is 



5Z. 



(0) 



5Z. 



(0) 



Zv = 0] ^ j^^jo) j 2 5. J (x, y) 6 J (x', y') 



(74) 



exp 



'(0) 



q (l + zfp^vrVA 



exp 



i + (2 + 4°Mp>V^ 



1 + 



1 + 4°V^vrVA^ 
The third term in Eq. ()28|) is 
1 6Z 



0s i^^y) <Ps i^'^y' 



V 



27(0) 



6'Z 



Zv[j = 0] 5ZP SJ{x,y) 6J{x',y' 



exp 



7(0) 



d 1 + 



vr 



27(0)" 



+ 



V 



exp 



0, (x', 2/0(1 + /) 



1 + 4°V^vrVA 
In k space the two contributions reads 

1 5^Zy bzf 6Z^'^ 



0, (x,y)0, (x', 2/0(1 + /) 



+ 



SZ^ 



V 



27(0) 



V [J = 0]s (z^^A ' -^-^ i^^ y) (^'' 2/0 Zv [3 = 0] ^40) 5 J (x, 2/) 5 J (x', y') 



5(k) 



1 + 4V^^v^ 



+ (1 + /) 1 + 



1 + 4V^vrVA^ 



/.^5(k) 2 + / + 



3 + / 



^/P^ + zP 



The results for Eq.(|28|) can now be collected 

1 , , 



vort+cs 



72 \V^FV^v)k 
^0 

' ^+ ^ 



(75) 



1 - 



p^P + 1 



k^pl \ k^plJ 87r2 \^ plk^ + 1 + zPAtx^pI/A^ 



+5(k) 2 + / + 



3 + / 



A/ Pi + ZP 
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We can make few remarks here. If the arbitrary position in plane of the 
vortices and the interaction between physical vortices were neglected, the 



only term that would persist is exp 



The first k-dependent factor in 



Eq.(??) comes from assuming that a statistical ensemble of realizations of 
the vorticity field is generated from the random positions in plane of the 
vortices, even reduced at a (5-type shape. In practical terms this may be 
represented as follows: in a plane, an ensemble of vortices can be placed at 
arbitrary positions. We construct the statistical ensemble of the realizations 
of this stochastic system. If we measure in one point the field, it will be 
zero for most of the realizations and it will be finite when it happens that 
a vortex is there. This is a random variable. Now, if we measure in two 
points and collect the results for all realizations, the statistical properties of 
this quantity (the two-point auto-correlation) has a Fourier transform that is 
given by the two factors multiplying the square bracket in Eq.(??), divided to 
(since we have the auto-correlation of the vorticity). When the interaction 
is considered, the factor in the curly bracket appears. 



5 Random field infiuenced by vortices with 
random positions 

Consider the equation 

F [0] = Vlifi - aip - f3ip'^ ^ (76) 
and extract from the total function the part that is due to the vortices 



N 



V {x, y) = ^ <t>^s^ {x, y) + (l) {x, y) 



(77) 



o=l 



Replacing in the equation we have 



N 



a=l 
N 



N 



-2/3 



a=l 



a=l 

(t> {x, y) 



N 



a=l 



(78) 



+Vl0 -a(j)- p(j)^ 
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The first line is zero and we have 



N 



a + 2/3j]0(")(x,y) 



a=l 



- ^ 



(79) 



We write a Lagrangean for the random field according to the MSR procedure 



>c[x,0] = x|vl0- 

and the action functional is 

S<pv [X, 0] = / dxdyx { 



N 



a + 2/3j]0(") {x,y) 



a=l 



(80) 



N 



a 



+ 2/3 ^ {x, y) 



a=l 



0-/302 



(81) 



y dxdy \ - (Vx) (V0) - X 



AT 



a + 2/3^0(") (a:,y) 



a=l 



The generating functional is defined from the functional integral 

j L'[x]i5[0]exp(-V[x>]) 



0-/3x0' 



(82) 



Now we will modify the action by considering as usual the interaction with 
external currents, 



5 [J] =M-'Jd [x] D [0] exp (- V [X, 0] + JxX + M) 



(83) 



With this functional integral we will have to calculate the free energy func- 
tional. The functional S [J] depends on the function representing the vor- 
tices. The vortices are assumed known but their position in plane is random 
therefore we have to average over them: —W [J] — (In (5 [J])). 

5.1 The average over the positions 

To perform the statistical average over the random positions of the vortices. 
The functional which we have to average is 



-1 



D[x]D 



(84) 



X exp jy dxdy [{Vx) (V0) + ax(t> + /3x0' 
+ (2/3f;0(«) (x,y)] x0 



a=l 
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The part that depends on the positions can be written 



exp I J dxdy \2(3 ^ {x, y)j x<P 

N 

2/?0^^x(ra)0(r„) 



(85) 



= ( exp 



a=l 



where (f)^ is now a simple amphtude. Consider the more general situation 
where we have to average in addition over the amplitudes 0* of the vortices. 
If 0^ is a stochastic variable we have to perform an average of the type 



exp 



N 



2/30^5^x(r„)0(r„) 

a=l 

= n(e^PWx(ra)0(ra)]>r„ 



(86) 



a=l 



Consider that the (now) random variable (ff has the probability density 

Then, restricting for the moment to only the average over 0*, 

(exp [2/30*x(ra)0(r„)])^, 

d<t>'g (0^) exp [2(34>\ (r„) (r„)] 



(87) 



-oo 
oo 



where 
Then 



d(t)'g {(f)') exp {iX(f)') 



iX (r„) = 2(3x (ra) (fa) 



(exp [2/30^x(ra)0(ra)]),. =5^[A (rj] 



(89) 



(90) 



where is the Fourier transform of the probability distribution function g. 
We use the notation 

g{vk)=g[X{va)] (91) 
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and we have to calculate 

N 



N 



a=l 0=1 

For this we introduce the function h 

g{ra) = h{ra) + l 
and we rewrite the average as 

N 



a=l 

N 

n(Mra) + l) 



(92) 



(93) 



(94) 



a=l 



J2 ^(rn)/i(rij.../i(rij 

;=0 ii<i2<...<ii 



exp 



°° 1 /■ 

il -I ^12 1 ■■■i^i 



where we have introduced the cumulants of the distribution of points rj^ in 
the plane. 

According to our assumption 



C(0 



1/A 1^1 
1>1 



(95) 



where A is the area in plane. Therefore we have the result of averaging 



(96) 



n(^(ra)) = exp ^ I drh{r) 

a=l L J . 

= exp 1^ J dr[g{r) - 1]| 

The part of the generating functional which depended on the positions can 
now be written 



exp 



N 



2/30^5^x(ra)0(ra) 



a=l 



expjiy dv[g{v)-l\^ (97) 
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and the full partition function is 



(98) 



X exp |y dxdy [( Vx) ( V0) + ax0 + /?X0^ 
+^ J dr[g{r)-l] | 



The action at the exponent is 



(99) 



The function g (r) has as argument the expression —i2(3x (r) (p (r). 

A reasonable assumption is that the vortices can only be positive or neg- 
ative and with the same magnitude. Then we have 



9 (0^) = o - 0o) + 5 if + 0o)] 



and the Fourier transform is 



9{q) 



df exp (^#^) - [6 (0^ -<Po) + S (0^ + 0o)] 



- [exp (ig0o) + exp (-ig0o)] 



(100) 



(101) 



= cos(g0o) 

This must be calculated for the argument —i2Px (r) (r) and gives 

^(r) = cos H2/30OX (r) (r)] 
= cosh[2/30ox(r)0(r)] 

The functional integral that must be calculated becomes 



(102) 



j D[x]D [0] 



(103) 



X exp 



dxdy 



( Vx) (V0) + q;x0 + /9X0' + -7 cosh (2/300X0) 



The system is perturbed with an external current J (x, y) acting on the field 

0(x,y). 



= J^-' j D [x] D [0] exp {Sj) 



(104) 
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/ 



dxdy 



( Vx) ( V0) + ax(l> + /3x0' + -r cosh (2/30oX0) + J0 



(105) 



6 Approximation for small amplitude vortices 



Consider that the amphtudes of the vortices are not high, 0q. Then 

cosh(2/30oX0)^l + ^(2/?0oX0)' 
and the action (removing some terms without significance) 



/ 



dxdy 



( Vx) (V(/>) + ax(l> + /5x0' + 



mo) 

A 



xV + ^0 



(106) 



(107) 



This action in principle can lead to a perturbative treatment but with 
two vertices, of order three and of order four, a very difficult and unusual 
problem. We remark however that 



dxdy 



-X ( + ax0 + PX<P^ + 



Wo) 



xV' + J<i> 



:io8i 



may become quadratic in x and in (f) (therefore the functional integral is 
Gaussian) if we succeed to separate the product x^0^- 



6.0.1 Technical step 

Consider the following formula to disentangle the two variables 



then 



exp ( ^2gx'0' 



and we have the action 
1 



dsexp ( — — •\/2Qx0* 



-1 



oo / ^ 

dsexT) I — 



X 



D [4>] D [x] exp \ / dxdy -x (V V) + ax0 + 



(109) 
(110) 

(111) 
2(5x0s + J(l> 
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or 



Sj^ J dxdy [J(l>] + J dxdy [-X + x (« - sy/2Q) + /3x0'] (112) 



We have used the notation 



Q = 



(113) 



The functional integration over x can be done immediately and gives a func- 
tional 6 with the argument the equation in the modified form. If we make 
the integration over x we obtain 



D [x] exp I J dxdy [-x (V^) + X - + Px<l>'] (|l4) 

= J D[x\exp J dxdyxF{(f)) 



.SF 

' 86 



(115) 



where we have introduced the notation 

F (0) = - V'(/> +(a- s^/2Q^ (j) + 

and the function (p^ is the solution of the differential equation F {(p) — i.e. 

F (^'^ = (116) 
Inserting this result in the integral over the functions 4>, we have 

.6F 



N'^ I D [0] exp ( / dxdyjA 



5 (117) 



.6F 



exp ( / dxdy J (f) 



Then 



1 



ex. / 



—I 



.SF 
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-1 



{/ 



exp ( / dxdy J (f) 



(118) 
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The most important part is to calculate 



SF 



(119) 



But this simplifies with the part of the factor 



(120) 



The fact that the initial factor is written before averaging over the positions 
and the second factor is obtained after averaging will make a small difference 
which can only be of higher order. What remains is a factor of normalization 
that we call c^^. 

The functional derivatives are 



6Zi 



1 1 



SJ {x,y) 
and 



c V27r 



J (is exp ^— -s^^ 0^ (x, y) exp dxdyJcj)^^ (121) 



1 1 



6Jix,y)5Jix',y') 



dsexp 



and taking 



c V27r 

x^"" (x,y)0^ {x\y')exp 
J = 



:i22i 



dxdyJcj)^^ 



(123) 



we can calculate the contribution to the two-point correlation arising from 
this part of the partition function. 
We have to solve 



+ (a - s^/2Q^ (f) + (3(f)'^ = 



(124) 



for a small amplitude field (p. This can be done by approximations, starting 
from the solution of the equation without random centers. 

We can see that the immediate effect of the scattering of the turbulent 
field by the vortices is a renormalization of the coefficient a of the equation. 

The solution that we need for Eq. ()124|) must refiect the fact that we 
are examining a turbulent field, in the dynamics of which the randomly 
placed vortices have been included. We can see that the integration over 
the auxiliary variable s must not necessarily be truncated to ensure that 
as = a — S\/2Q remains positive. This is because the field is turbulent and 
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there is no restriction concerning the existence of the vortex solution which 
should be imposed to the turbulent component of the field. We also note 
that the approximation taking only the s = contribution in the integral 
()122j) returns us to the problem of the turbulent field without interaction 
with the random vortices. Although the main contribution (as shown by the 
Gaussian integration) comes from the free turbulent field itself, the presence 
of vortices is contained in the rest of the integral, for |s| not close to 0. 

These remarks suggest to study the statistical properties of the trans- 
formed field, where a is replaced with a^. 

7 The background turbulence: perturbative 
treatment 

The vortices represent the strongly nonlinear part of the system and they 
have been extracted and treated separately as shown in Section IV. How- 
ever, there is still nonlinear interaction in the remaining random field. This is 
the nonlinear mode coupling which produces the stationary turbulent states 
of the system, even in absence of any definite structure formation. This 
interaction must be taken into account when we analyse the statistical prop- 
erties of the turbulent field. From general consideration we know what can 
be expected from this analysis: we will obtain a small departure from pure 
Gaussian statistics, expressed in nonlinear renormalization of the propagator 
and of the vertex of the interaction. The presence of random interaction with 
elements of the system that are beyond our simple model is accounted, as 
usual, by a noise term acting like a drive for the system of random waves. 
The noise will be assumed with the simplest (white) statistics 



and should be considered as a random stirring force composed, for example, 
of random growths and decays of marginally stable modes, thus injecting at 
random places some energy into the system. 

We will provide a standard perturbative treatment for the differential 
equation 



with the objective to calculate correlation functions and other statistical 
properties. 

As usual we start by defining the generating functional of the statistical 



(125) 




(126) 
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correlation 

exp{-W[J]) (127) 
= Jv[r)]V [(f)] exp I J dxdy [-77 + 77q;,0 + 77/30^ + J^Ty + J^0] 

x(^Jv[r)]V [0] exp | ^ dxdyTyC 



where we have introduced the notation 

as = a- s^2Q (128) 
The second factor can easily be calculated 

jv[ri\V [0] exp j dxdyriC^ ^ (129) 



V [t]] V [(j)] exp 



dxdy {2Drf) 



In the generating functional there are the two dual functions, r] (x, y) and 
(j){x,y), the second being the physical field of the random waves. We have 
inserted the action related with interaction with "external" currents, J,, and 
Jfj) which will permit to express the averages as functional derivatives. 
Defining the Lagrangean density (after an integration by parts) 

jr ^ (Vr/) (V0) + rias(p + vP4>'^ + Drf" + J^ry + J^cf) (130) 

we will separate into the Gaussian the nonlinear interaction part 

C^Co + Ci (131) 

£0 = (Vt;) (V0) + r]as(l> + Dri" + J^ry + (132) 

Ci = rip(l)^ (133) 

The first part of the Lagrangean is not linear because of the square term, 
but it is Gaussian therefore it would not pose any problem to the functional 
integration. We note that we have an order three vertex. 

The action functional is correspondingly divided into two parts and can 
be written, as usual 

Z [J] = Jv [r]] V [0] exp (^J dxdyC^ exp {^j dxdyC^ (134) 
= exp j dxdy SJ^ Jll) j M ^ exp dxdyCo 
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For the linear part we have the Euler-Lagrange equations 

d 6Co ^ d 6Co _ 6Co ^ ^ ^^^^^ 



dx S {drj I dx) dy S {drj /dy) Srj 

d SjCq d 5jCo SjCo „ 

+ — ^ /o / /o \ FT ^ ^ [^'^^) 



dx5{d(l)/dx) dy5{d(f)/dy) Scf) 
The equations are 

- - 2Dri = (137) 

and 

A77 - a,r] = (138) 
The latter is an inhomogeneous Helmholtz equation and has the solution 

Vo {x, ^) = / dx'dy'G^^ (x, y; x\ y') (x', y') (139) 

in terms of the Green function appropriate for the space domain of our anal- 
ysis and taking into account the boundary conditions for 77. 
The first Euler-Lagrange equation gives 

A0-a^0 = J^ + 2Dr]Q (140) 

= Jr, + 2D J dx'dy'Gr,4,{x,y;x',y') J^{x',y') 

with the solution 

00 = J dx'dy'G4,r,{x,y;x',y') Jr,{x',y') (141) 

+2D J dx'dy'dx"dy"G^,{x,y;x",y") G,^ {x",y";x',y') J^{x',y') 
We can introduce the matrix of propagators 

G = (G,,) = ^^^^ (142) 

where 

G^^ {x, y; x', y') = 2D j dx"dy"G^, {x, y; x", y") G,^ (a;", y"- x' , y') (143) 
Introducing the notation 



; ) (144) 
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we will write symbolically (summation over repeated indices is implied) 

^. = I,.WG..J. (145) 



for 



dx'dy'G {x, y; x', y') Jk (x', y' 



(146) 



We can now calculate the linear part of the action along the system 's 
configurations given by these solutions. The action will be, of course, a 
functional of the current. 



Soj = j dxdy [(Vr^o) (^^o) + «s^o0o + ^vl + Jr,Vo + J'P'Po] (147) 
= ^ y dxdy [-r^oA^o + a.r/o^o + ^Dril + J^ri^ + J^cpo] 

+^ J dxdy[-(j)QAr]Q + asr]Q(f)o + Jr,Vo + J<f<l>o] 



or 



Soj = ^ j dxdy [rjo (-A0o + "s^o + '^^Vo + Jr^) + J,f'Po\ 
+ dxdy [00 (-Ar/o + + J^) + -^^y^] 

We now express the solution by the Green functions 



Soj ^ 2 ' ^^^^ 



j dx'dy'G^ {x, y; x', y') {x' , y') 

+Jr, J dx'dy'Gr,4, {x,y;x',y') J^{x',y') 
which can be written using the convention introduced above 



We return to the generating functional 

S S S 



Z[J] 



exp I /3 / dxdy 



exp ( / dxdyC, 



5Jrf 5J^ 5J^ 
5 5 5 



exp {Sqj) 



ijk 



5Ji 5Jj 5Jk 



X exp 



(148) 



(149) 



(150) 



(151) 



^ J dxdy J dx'dy'Ji{x,y)Gij{x,y;x',y') Jk{x',y') 
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where we have introduced, for uniformity of notation 



Cijk {x, y) = PSirjSj^Sk^ (152) 

This is a classical framework for diagrammatic expansion. 

In general the functional approach is well suited for statistical problems 
where the statistical ensemble of the realization of the system's configurations 
is produced by the effect of a noise or by a random choice of initial condition, 
as explained in MSR and Jensen. In the most usual case, where there is a 
noise acting on the system, using the (known) statistical properties of the 
noise means implicitly that we make a change of variables, from the field 
to the noise C and this introduces a Jacobian, 

J[<l>]=V[C]/V[<j>] (153) 

The Jacobian cancels all the diagrams that correspond to the vacuum to 
vacuum transitions, or equal-coordinates correlations. For this reason it is 
not mentioned. This is explained in Appendix B. 

When applying the operator to the second exponential we have to remind 
that we need at least two open ends with currents since we intend to 
determine the spectrum. On the other hand, it is clear that any open end 
with isolated either or would vanish after taking the currents to zero. 

The first term is naturally the diffusion driven by the random rise and 
decay of marginally stable modes, represented here as a noise. 

The next term that is useful consists of the product of two vertex-like 
operators applied on a product of four propagator-like terms 

Obviously, this is the one-loop diagram. We will derivate to the first and last 
factors (currents) so that, more explicitly, its structure is 

G(j)rfC^(l)(j,G(j)(j)G(i,(i,C(i,(jy,^G^fl) (155) 

We can represent the Green function by its Fourier transform 

(A - as) G (x, y- x\ y') = -5 {x - x') 5 {y - y') (156) 



G{x,y;x',y') = J 



dkxdky-^ 



1.2 



(157) 



X exp [-ikx {x - x')] exp [-iky {y - y')] 
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We obtain the propagators for the couphng between the field 4> and its dual 
G^r,{x,y;x',y') = [ dk^dky }^ — (158) 

J K,^-\- K,y -\- as 

X exp [-ik^ {x - x')] exp [-iky {y - y')] 
= Gr^^(x,y;x',y') 

G^^{x,y) = 2D J dx"dy"G^,{x,y;x",y")Gr,^{x",y";x',y') (159) 



"{kl + kl + asY 
X exp [-ik^ {x - x')] exp [-iky {y - y')] 

The intermediate integration over space {x", y") yields equality of the Fourier 
variables of the two Green functions. We then have 

G^^ = 2D {kl + kl + as)~^ (160) 

The one-loop term is expressed as two integrals over the intermediate 
momenta at the two vertices where we have a product of three Green func- 
tions with the vertex. Conservation of momentum in the loop and overall 
conservation of the diagram (which means that the k at input line must be 
the same at the output line) lead to 

(<t,(x,y)<t,(x',y')) (161) 

dM*, W), exp 01 
2D 



(k2 + asf 

[ [ I 2D 2D 1 

+ dkV dp— (3- -2 2/3- 



k^ + a/ (p2 + asf[{k-pf + asf^' + 

The coefficient F represents the multiplicity of this diagram and factors of 
normalization. 

(</>0)k (162) 

\2 



^i^Dpy TIT— dp 



M' + aJ J "(p2 + a,)' [(k-p)2 + a,] 
= r(2i?/3)^(-— ) I dp 



1 V f , 1 1 



+ y J [{k/2 -pf + a^Y [(k/2 + pf + asY 
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We transform the integral 



1 



[p2-k-p + kV4 + Q;,]- 
1 



[p2 + k-p + kV4 + a,]' 
Now we introduce the variables that will replace the constant terms 

b = k74 + a, 

and 

k ■ p = kp cos 9 

= up 

The integral becomes 

1=1 pdpdO 



\p'^ — up -\- 6] [p2 + up + V\ 
We write the denominator in the form 

[p^ — up + h\^ [p^ + lip + 6] ^ 

= \p^ + hf-uY\ 

= [p^ + 2p% + h''-uYf 

= [p'+{2b-u^)p' + h^Y 
= (p^ + 27p2 + (5)' 



where 



27 = 2h-u^ 
5 = 



The integral becomes 



pdp 



(p4 + 27p2 + 6y 



we make the substitution 



p ^ X 
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1 /-oo 

de- 



dx 



2 Jo (x2 + 2-fx + 5f 
The integration over the intermediate momentum is 

d.v 



Here 



{x^ + 27a; + 5y 



(171) 



(172) 



2 {5 - ^^f'^ 



TT 7 

— — arctan — , 
2 



+ 



V 2(5/ (5 - 72 



7 = 



= - [A;^ (1 -2cos^e) +4q;,] 



(173) 



7 



6 - (/72 



1 A:2(l -2cos2g) +4Qi, 

2 A; |cos6'| a/ A;2 sin^ 6* + 4q;s 



(174) 



Then 



= cos^ ^ sin^ ^ + a,) 
7 _^A;2(1-2cos2^)+4q;s 



(F + 4q;,)^ 



pdpd9- 



(175) 



[p2 - k • p + k74 + a,]' [p2 + k • p + kV4 + a,] 

- /-{ 



(176) 



.(177) 



X 



cos2^(A;2sin2^ + Q;,)]^^^ 

TT fc2(l -2cos2^) +4a, 
9;rc1j9;n ^ 

2 2A;|cos^| sin2^ + 4a. 



fc2(l-2cos2^) +4a. 



A;2cos2^(fc2sin2^ + a,) (A;2 + 4a,)^ 
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We can show that the integral is not singular at the limit 



cos6' 







(178) 



so that the integration over 6 is safe. 

As order of magnitude the integral is dominated by terms like 



while the diffusive part i 




(179) 



2D 



(180) 



After this we obtain for the turbulent background field 




(181) 



8 Summary 

It is usually assumed that a turbulent plasma (in prticular with embedded 
structures) should exhibit a spectrum of exponential or algebraic type. There 
is no universal theoretical basis for such an assumption except for cases where 
the scaling invariance is justified on physical ground. We find that it would 
be more adequate to extract a particular behavior (on different spectral in- 
tervals) from expressions like Eq. ()75|1 and Eq. ()181|) , via regression on simple 
numerical values. This can be done for a particular physical system. Here, 
however, we will return to the traditional exponents in order to exhibit some 
associations we consider to be general. 

We can now collect the results of the analysis. In the left column we write 
the approximative behavior of the contributions and in the right column the 
physical origin, according to this theory. 



k-' 



~ S (k) k-^ (weak pI/A) 

S (k) [1 + / (k)] 
/,^onm2 1 1 
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We see two ways in which these results can be useful. 

First, the combination of various exponential dependences from the listed 
contributions will result in an overall dependence with exponents that may 
be compared with the experiments or numerical simulations fHl, [201, [21], 
|22j . |24j . The weight of each contribution is dependent on factors like the 
shape of the vortices (including the spatial extension compared to the area 
A), amplitude of the vortices, strength of the random drive (D). Also in 
the problem enters as parameter the function of distribution on the plane 
(assumed here uniform). The energy of interaction between vortices may be 
reconsidered along the example of vortices in superfluids. 

Second, the spectrum can be seen as dominated in different spectral do- 
mains by one or another of these contributions. This is compatible with the 
known fact that the spectrum has regions with different exponents /i. 

In both these ways we have to solve an inverse problem, starting from the 
experimental (or numerical) spectrum and using the above list to map the 
exponential form to a physical process which may have been at its origin. 

Various extensions of this treatment are possible. One can introduce a 
chemical potential for the description of the statistical equilibrium consisting 
of generation and suppression of coherent vortices. There are treatments of 
this type for the conversion of the global rotation of superfluids into localised 
vortices, an example that may also be useful for the consideration of the zonal 
flow saturation in tokamak, besides the Kelvin- Helmholtz instability and the 
collisions. 

The treatment by generating functional allows in principle determination 
of statistical correlations at any high order desired. However, while the func- 
tional derivatives can easily produce the n-th order cumulant, we must be 
sure that the generating functional has been calculated with the necessary 
precision, or, in other words, we must be sure that we have incorporated the 
physical origin of these correlations. This method should be accompagned 
by the more standard analysis of closure of the hierarchy of equations for the 
correlations and these two approaches must be seen as complementary. 

As a final remark, the physical model adopted in the present treatment 
may be extended to cover more complex regimes, with, of course, a certain 
increase in the analytical work. 
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9 Appendix A : physics of the equation 



It is of interest to study a realistic two-dimensional model, like Hasegawa- 
Wakatani or similar. However we need for the beginning a simpler equation 
and if possible with a vortex solution with known analytical expression. 

A possibility is the equation for the ion drift instabilities but here we 
must specify the scales. 

Consider the equations for the ITG model in two- dimensions with adia- 
batic electrons: 

^ + V-(v,n,) = (182) 

(9vj , , e , ,x e ^ 

-l + {Y,.V)vi = — -V0 + — ViXB 
at rrii rrii 

We assume the quasineutrality 

Hi a; Ue (183) 
and the Boltzmann distribution of the electrons along the magnetic field line 

ne = noexp|^-i^^ (184) 
In general the electron temperature can be a function of the radial variable 

Te = Te {X) (185) 

The velocity of the ion fluid is perpendicular on the magnetic field and is 
composed of the diamagnetic, electric and polarization drift terms 

= vxi (186) 

— ^dia,i + + ^pol,i 

Ti 1 drii^ 



\e\ B rii dr ^ 
— V(/) X n 



+ 

+ (v£ • Vx) Vx0 



BQi \dt 

The diamagnetic velocity will be neglected. Introducing this velocity into the 
continuity equation, one obtains an equation for the electrostatic potential 
0. 
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Before writing this equation we introduce new dimensional units for the 
variables. 



iphys 



e 



T 

e 



phys 

(187) 



r^phys yphys 
Ps ' Ps 



tphys ^ ^phys^^ ^^gg^ 

The new variables {t, x, y) and the function are non-dimensional. In the 
following the primes are not written. 

9.1 The equation 

With these variables the equation obtained is 

I (1 - Vl) (190) 
- {-V±(f) xn) -Vd 
+ (— X n) • 
+ [(- Vx0 X n) ■ V J (- Vi0) 
= 

where 

Vd = -V^lnno- V^lnTe (191) 
vr = -V_LlnTe 

(This is Eq.(8) from the paper [TU]). 

9.2 No temperature gradient 

From the various versions of the nonlinear drift equation, (in particular 
the Hasegawa-Mima equation ) we choose the radially symmetric Flierl- 
Petviashvili soliton equation: 

(1-.M)|..|-..|^0 (192) 

where = Cg/^i, Cg = {Te/mi)^^'^ and the potential is scaled as = ^ . 
Here L„ and Lt are respectively the gradient lengths of the density and tem- 
perature. The velocity is the diamagnetic velocity Vd = ^j^- The condition 
for the validity of this equation are: (kxPg) [kp^Y ^ 'He'iT^ where rj^ = j^. 
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The exact solution of the equation is 
1 ) sec h"^ — 



{y,'t;yo,u) 



1 



Vd' 

u 



\ 1/2 

j {y-yo- ut) 



(193) 



where the velocity is restricted to the intervals u > Vd or m < 0. In the 
Ref. jHj the radial extension of the solution is estimated as: (Ax)^ ~ p^Ln. 
In our work we shall assume that u is very close to , m > (i.e. the 
solitons have small amplitudes). 



10 Appendix B : Connection between the MSR 
formalism and Onsager-Machlup 

In our approach the most natural way of proceeding with a stochastic differ- 
ential equation is to use the MSR type reasoning in the Jensen formulation. 
The equation is discretized in space and time and selected with 6 functions 
in an ensemble of functions (actually in sets of arbitrary numbers at every 
point of discretization). The result is a functional integral. There is however 
a particular aspect that needs careful analysis, as mentioned in the previous 
Subsection. It is the problem of the Jacobian associated with the 5 functions. 
This problem is discussed in Ref. i34j and for consistency we include here the 
essential of their original treatment. 

The equation they analyse is in the time domain and is presented in most 
general form as 

where the number of stochastic equations is , if is functional of the fields, 
Vj is the streaming term which obeys a current-conserving type relation 

^\/,Mexp{-ifM} = 

The noise is 9j. 

The following generating functional can be written 

Ze = j D [0,. it)] exp j dt {t)] \[ 6 (^-^ + K, [0 (t)] - 9^ J [0] 



the functions Ij (t) are currents. 
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and J [0] is the Jacobian associated to the Dirac 5 functions in each point of 
discretization. 

The Jacobian can be written 



J = det 

Up to a multiphcative constant 



d_ 5K^ 



or 



J = exp Tr In 



J = exp I Tr In 



d_ 5K_ 



S{t- t') 



6 {t - f 



dt 



5{t- t') 



1 + 



5K jt) 



Since the operator (^) ^ is retarded, only the lowest order term survives 
after taking the trace 



J = exp 



dt 



SK, [0 jt)] 
54^ At) 



The factor 1/2 comes from value of the © function at zero. 

In the treatment which preserves the dual function associated to in 
the functional, there is a part of the action 

4>K[<t>] 

Then a and a of the same coupling term from (f)K [0] close onto a loop. 
Since G-^^ is retarded, all these contributions vanish except the one with a 
single propagator line. This cancels exactly, in all orders, the part coming 
from the Jacobian. 

Then it is used to ignore all such loops and together with the Jacobian. 

In conclusion we can compared the two starting points in a functional 
approach: The one that uses dual functions (t) and x (t), closer in spirit to 
MSR; And the approach based on Onsager-Machlup functional, traditionally 
employed for the determination of the probabilities. Either we keep x (t) 
ignore the Jacobian (the first approach) or integrate from the beginning over 
X (t) and include the Jacobian. The approaches are equivalent but we have 
followed the first one. 
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